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A procedure is described for matching a given stationary axisym- 
metric perfect fluid solution to a not necessarily asymptotically flat 
vacuum exterior. Using data on the zero pressure surface, the proce- 
q-i dure yields the Ernst potential of the matching vacuum metric on the 

5^ | symmetry axis. From this the full metric can be constructed using a 

variety of well established procedures. 

> 

u • 1 Introduction 

Because of the tremendous strides that were taken during the 1970's and 
1980's in coping with the mathematical problems presented by stationary 
axisymmetric vacuum (SAV) and electrovac fields, the attention of many 
workers in this field has shifted in recent years to those even more challeng- 
ing problems that are associated with stationary axisymmetric perfect fluid 
(SAPF) spacetimes. Except in the case of zero pressure fluids (dust) there 
exists no general method for solving the Einstein field equations within the 
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fluid itself. Neither has any general method for joining a given SAPF solu- 
tion to a vacuum SAV exterior at a zero pressure surface been described. It 
is the latter problem to which we shall propose a solution in this paper. 

From the outset we stress that we are not insisting that the SAV metric 
be asymptotically flat. Ours is a method with which to effect a matching 
of a given SAPF to a SAV exterior, and not a method for solving a global 
problem. In the latter case, the asymptotic flatness and the singularity free 
nature of the solution are built in ab initio, and the SAPF metric as well as 
the SAV metric are determined a posteriori. The rotating dust disk solution 
of Neugebauer and Meinel|I[] remains the only known global solution of the 
stationary axisymmetric field equations corresponding to a physically inter- 
esting rotating source, the analog of the classical zero pressure MacLaurin 
disk. 

Restricting their attention to dust disks, Neugebauer and Meinel em- 
ployed a method involving two integral equations, the "small" and "big" 
integral equations. The former integral equation was derived from the Neuge- 
bauer linear system for the Ernst equation under the assumption that on the 
disk < p < p the Ernst potential £(0,p) is real and independent of p, 
and under the assumption that there exists a singularity free asymptotically 
flat matching SAV metric. The solution of this first integral equation yielded 
the axis values S(z, 0) for z > of the Ernst potential. The second integral 
equation was derived from a Riemann-Hilbert problem, and was used in or- 
der to obtain the Ernst potential S(z, p) of the vacuum metric that matched 
the axis data derived in the first step. Finally, the surface density of the dust 
disk was determined by computing the jump in one of the metric fields as 
the disk is crossed. 

In the present paper, and with an eye toward providing a procedure that 
can be employed routinely by those who discover new perfect fluid solutions, 
we develop the basic machinery that will be required to treat the "joining 
problem" for a broad class of spinning fluid spheroids that will be defined in 
Sec. [3p. Our approach, which is a variation of the "monodromy transform" 
approach of Alekseev0, employs matching data, essentially certain metric 
components and their differentials, on the zero pressure surface to compute 
a 2 x 2 matrix function II (r) of a complex spectral parameter r, the so-called 
"monodromy matrix." 

The axis values of the Ernst potential of the exterior SAV solution can be 
determined immediately from II (r), which means that any of the standard 
methods of generating a SAV solution from its axis data can be employed. 
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Thus, for example, the SAV solution can be constructed by solving a homo- 
geneous Hilbert problem (HHP) that was developed by the authors]^, [|. HI 
to effect Kinnersley-ChitreH transformations. 

Since we are dealing with a group, there exist a number of sequences of 
transformations that can be used in principle to transform Minkowski space 
into the final SAV metric. Certain transformations may be found to be eas- 
ier to carry out than others, so, here, as in many things, experience is the 
best teacher. However, if one wants a standardized treatment method, that 
would be obtained by employing an HKX transformation 0, either carried 
out as described by the authors of the HKX transformation, or by transform- 
ing Minkowski space into a static Weyl metric, a step that requires only a 
quadrature, then transforming that Weyl metric into the final SAV, a step 
that involves solving an ordinary Fredholm equation of the second kind[||. 
Moreover, a possible alternative approach has been provided by Klein and 
Richter||. 

Even if the n-matrix cannot be evaluated in closed form, by studying its 
domain of holomorphy in the extended r-plane one can ascertain whether or 
not an asymptotically flat singularity free global solution exists. In Sec. |] 
we shall illustrate the various steps of the n-matrix method with two exam- 
ples. The Schwarzschild interior solution is ridiculously simple but has the 
advantage of involving only short and easy calculations. The second exam- 
ple involves a rigidly rotating stationary axisymmetric dust spacetime with a 
dust density distribution that is specified a priori. There, since the pressure 
vanishes identically, one has the option of selecting at will the matching sur- 
face. After choosing a simple surface, we were able to show that the matching 
SAV metric is certainly not asymptotically flat. It would be exciting if one 
were to find a dust density distribution and a choice of matching surface that 
resulted in an asymptotically flat matching SAV metric. 

Constructing a matching vacuum metric for the Wahlquist interior solu- 
tion is likely to be challenging. In that case we do not expect the matching 
SAV spacetime to be asymptotically flat. We hope that the availability of 
the general procedure for matching a given fluid solution to a vacuum met- 
ric that we shall describe in this paper will encourage those who discover a 
new exact SAPF solution to investigate possible physical interpretations of 
their fluid solution, even if it turns out that the matching SAV metric is not 
asymptotically flat. 
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2 The Perfect Fluid Field Equations 



In this section we shall formulate the field equations that will be assumed 
to hold both within the region occupied by a perfect fluid and within the 
vacuum region outside the source, and we shall stipulate our continuity- 
differentiability premises. 
Let 

hcd :=K C -K d (c,d€ {3,4}), (2.1) 

where K 3 and K 4 denote the usual rotational and timelike Killing vector 
fields, respectively, whereupon h 33 > except on the axis, and /i 44 < 
everywhere. Let (Ai 2 ,g^) denote any one of the usual two dimensional 
Riemannian subspacesQ of the spacetime such that its Hausdorff space Ai 2 is 
orthogonal to the Killing vector orbits, and its atlas is C°°. All the differential 
forms that we employ in this paper will be on a two dimensional Hausdorff 
space Ai 2 or on subspaces thereof. In particular, /, u> and P will denote the 
0-forms with domain Ai 2 and values 

/(x) = -Mx), "(x) := Mx)//(x) (2.2) 

and 

P(x) := [/i 34 (x) 2 - /i 3 3(x)Mx)] 1/2 > (2.3) 

for all x G Ai 2 . Furthermore, eu will denote our choice of the unit 2-form 
on Ai 2 , and * will denote the duality operator for the differential forms on 
Ai 2 such that 

*e i2 = 1, *1 = e 12 (2.4) 
and, for all 1-forms v and w on Ai 2 , 

(-kv)w = e 12 (v ■ w), (2.5) 

where 1 denotes the identity mapping on At 2 , and we follow the practice of 
suppressing 'A' in all exterior products and derivatives of differential forms. 
One readily shows that, for all 1-forms v and w on Ai 2 , 

■k-kV — —V, v(kw) = w(kv). (2.6) 

Note thartfl 

P(x) = if and only if x £ M 1 ^, (2.7) 

where Ai\ x denotes the axis subspace of Ai 2 . 

1 These are the surfaces of constant (f,t). 

2 This excludes event horizons from consideration. 
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A The Matching Premises 

In this first exposition of the Il-matrix formalism, we shall restrict attention 
to those Ai 2 , the topology of which is described by Fig. 1.0 We assume that 
the common boundary Ai^ of the fluid interior Ai 2 in and vacuum exterior 
Ai 2 x has a C 1 parametric equation 

x = x(A) (0< A<tt) (2.8) 

that defines a homeomorphism of [0, ir] onto Ail such that the "north" and 
"south" poles, where the zero pressure surface Ail an d the axis Ai\ x inter- 
sect, are given by 

n = x(0), s = x(tt). (2.9) 
We also adopt the following (often tacitly assumed) premises: 

For all C n 1-forms v and w on Ai 2 , *v and v ■ w are C 1 if n > 1 
and C° if n = 0. Also, p is C°, / and u are C 1 , while P is C 2 , . 
and dP(x) ^ at all x G Ai 2 . Finally, there exists a continuous 
0-form A such that PA is C 1 and u = P 2 A. 

These premises include all of the matching conditions at the zero pressure 
surface A^q, and they are consistent with the field equations that we shall 
now specify. 

B Four Field Equations and the Euler Equation 

It is convenient to introduce the 0-forms 

r := -~\n(dP ■ dP) and ^ := f (2.11) 

and the 1-forms 

f := -P- l f 2 * duj := -f * [AdP + d(PA)] (2.12) 

and 

r] := P-kdnj) - -u£. (2.13) 
3 This restriction rules out, for example, sources of a toroidal shape. 
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In view of our premises, T and ip are C 1 , and £ and 77 are C°. From the 
definition (|2Tl| ) of T, 



e 1 :=e T i<dP and e 2 := e r dP (2.14) 
constitute an orthonormal pair of 1-forms on jVl 2 , and 

e 12 = ei e 2 = e 2V (*dP)dP. (2.15) 
We shall also employ the 0-forms £/, g?/0 an d ^/i"^ f° r which 

£ = * dP + (£//) and dip = (d^) * dP + (d n ip) dP. (2.16) 

The restrictions to the fluid interior M. 2 n and the vacuum exterior M. 2 X 
of each of the differential forms on M 2 that we have defined above will 
be denoted by the letter employed for the differential form with the affixed 
subscripts 'in' and 'ex', respectively. 

We now introduce a function a that will be called the boost form on M. 2 n . 
Let U denote the fluid world velocity field on Aif n , where we recall that 
U • U = — 1. Another timelike unit vector on M.\ n is the restriction to M.\ n 
of 

e t := -^j K 4- (2-17) 

It is easy to show that \e t ■ U| > 1 throughout M.\ n . We select the sign of 
K 4 so that e t • U > 1 and let a denote the 0-form on M 2 n such that 

cosha(x) := e t ■ U(x), sgn a := sgn U 3 (2-18) 

for all x G M. 2 n . The boost form will be employed in our expressions for the 
field equations. 

To ensure that the field equations are meaningful at all points of M. 2 n U 
M. 2 ex including points on M 1 ^, we adopt the following premises, which are 
not necessarily independent of one another: 

For all C°° 1-forms v and w on A4 2 n U A4 2 x , *v and v ■ w are C°°. 
Also, fi n , f ex , Pi n , P ex , A.^, A-exi Piny ^in an d P in ^ sinh a are C°°. 

(2.19) 

These premises are certainly stronger than necessary, but they enable us to 
avoid complications that would obscure this first exposition of the Il-matrix 
method. 
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With the above premises, all but one of the field equations are as follows: 

d(T + V) -P {(*dif))dri/) + (dil>)d n i/; 

+(2/)- 2 [(*£)£j +(£)£//]} = -KP P e 2F dP, (2.20) 

d*dP = -2nPpe 12 , (2.21) 
d£ = K f(p + e)(sinh 2a) e 12 , (2.22) 

d V = --nP[2p+{p + e)cosh2a}e 12 , (2.23) 

where k := 8nG/c 4 in conventional metric units, and it is to be understood 
that the restriction of the right sides of Eqs. ( p. 22] ) and ( |2.23|) to M. 2 X are 
identically zero. The contracted Bianchi identity in M.\ x is identically sat- 
isfied, and in A4f n it is equivalent to the following general relativistic Euler 
equation (in which we suppress the subscripts 'in'): 

dp+(p+e)[(cosh2a)#-(cosh2a-l)dP/2P-(siiih2a)*^/2/] = 0. (2.24) 



Note that the field equations ( |2.20|) and (|2.21|) hold throughout M 2 in 



eluding M l . From Eqs. flTfj) and ( |2T20| ), T(x) + ^(x) = T(n) + ^(n) at all 



x G M l ax . As is well known, the condition 

T(x) + ^(x) = at all x G M l ax (2.25) 

isQ necessary and sufficient for g to be locally Minkowskian on Ai 2 x ; i.e., 
to have no "conical singularities" on M. 2 ax . Therefore, we shall henceforth 
assume that r(n) + ^(n) = 0. 

Unlike Eqs. O) and (|2~2ll) , the field equations (^22|) and are 
generally defined only in the domain A4f n U M. 2 ex . This situation changes 
when e(x) = at all x G M.\, but we shall not pursue that topic in this 
paper. 

It is desirable to extend the domains of all of the field equations in M? n to 
M 2 n := M 2 n UMl, of all of the field equations in M 2 ex to M 2 ex := M 2 ex UMl, 
and of the Euler equation ( [2.24j ) to M.f n - To accomplish this, it is sufficient^ 



4 granted our premises 

5 It is to be understood that all derivatives at boundary points of M\ n and of }A\ :L 
defined using only sequences of points in M\ n and M-1 X , respectively. 
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to introduce the following reasonable premises, which are consistent with all 
other premises of this paper and with Eqs. (|2.20|) to (|2.24 ): 



1-forms v and w on A4 2 n , 



For all C 2 

are C , and Pi nj fini A{ n P{ n ^ Aj n , 

Pr 1 sinh a have C 3 , C 2 , C 2 , C 1 , 
extensions, respectively, to M. 2 n - 
v and w on M. 2 , *v and v ■ w are 



•kv and v ■ w 

Pini tin an d 

C 1 , C° and C° 
For all C°° 1-forms 

C°°; and P eX} f ex 



and A PT have C°° extensions to M. 2 



(2.26) 



The only remaining field equation (one which we have chosen to suppress 
in this paper) happens to be a 2-form equation involving d-kdT. A well known 
proposition asserts that, if (f ex , P ex , A ex ) is a solution of Eqs. ( |2.21| ), ( |2.22| ) 
and Q2.23|) on M. 2 ex1 then 6 Eq. ( [2.20|) on M. 2 X is completely integrable; and, if 
T ex is the integral of Eq. Q2.2Q ) on M 2 exL then (f ex , P ex , A ex , T ex ) identically 
satisfies the d * dY field equation on M. 2 X . It can also be shown that if 
{f in ,P in ,A in ,T in ,p im e in ,a) is a solution of Eqs. (CT , ( CT , , CT ) 
and (suppressing the subscript 'in') 



di(r + V) - 2P {c^ dj,^ + (2/)- 2 £/ Cu} = o 



on .M 2 n , thenF] Eq. ( |2.20|) on M. 2 n is completely integrable 



the integral of Eq. fl2.20|) on M^ n , then (/ in , Pj 
satisfies the d * c?r field equation on A4 



i at Ai n) Tj n , Pin, 



(2.27) 

and, if T in is 
:j n ) identically 



t 2 . 



C The Vacuum Region 



We now focus attention on Eqs. ( ggg ), ( ggg ) and ( gjg ) on .A/1 2 X . Since A^ 2 X 



is simply connected, these field equations are equivalent to the statement^ 
that C°° 0-forms Z, x an d exist on M. 2 X such that 

*dP = dZ, (2.28) 



6 Note that we are taking the liberty of employing the same notations for differential 
forms on -Mf n and M.\ x as we do for their restrictions to .Mf n and A4g X , respectively. 
However, keep in mind that the differential forms that we are considering are C°° in 
Mf n \J M 2 ex , but are not necessarily C°° on Mq. 

7 Ibid 

8 We are taking the liberty of suppressing the subscripts 'ex' in 'P ex \ l £,ex' and 'r) ex ' . 
Similar abbreviating liberties will be employed in later equations, and we shall depend on 
the context to help the reader avoid confusion. 
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V 



dx, 



(2.29) 
(2.30) 



Equation ( |2.28| ) enables us to introduce the Weyl canonical chart x 
(z,p) : = (Z(x),P(x)) that maps M. 2 ex onto 



{(Z(x),P(x)):xe^J- 



(2.31) 



As is well known, Eqs. ( |2.29| ) and ( p.30| ) then yield the following elliptic Ernst 
equation for the potential £ := f + %% expressed as a function S(z,p) = 
f(z, p) + ix{z, p) of Weyl canonical coordinates: 



+ P 



d£(z,p) 
dz 

d£(z,p) 
dz 



d_ 

dp 



+ 



d£(z,p) 
P dp 

d£(z,p) 
dp 



(2.32) 



throughout T> ex . We shall describe a way to determine the axis values of 
this ^-potential from the matching data at the zero pressure surface. Many 
methods are known that permit the construction of £(z,p) from S(z,0), 
assuming that suitable premises are satisfied. It is, however, not our purpose 
to rehash these well known procedures. 



3 The Monodromy Matrix II (r 

A Determination of the Matching Data 

The curve 

P :=p(x),P(x)):x6M;} (3.1) 
that represents the matching surface is given by the parametric equations 



z(\) := Z(x(A)), p(X) := P(x(A)) (0 < A < vr), 

and the matching data comprise z(X), p(A) and 

/(A) := /(x(A)), xW ■= x(x(A)), 
w(A):=w(x(A)), 0(A) := 0(x(A)). 



(3.2) 



(3.3) 
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The functions of A that are defined above are to be determined from the 
given fluid solution in fA\ n as follows. 

The given fluid solution is a compact two-dimensional Riemannian space 
(Ai 2 ,g^ ) with a maximal C°° atlas. If (x 1 ,^ 2 ) denotes the coordinate pair 
corresponding to any point x in the domain of a chart in this atlas, one can 
use the metric gf^ to compute the 0-forms for which 

*dx* = dx J y)(*). (3.4) 

Then we may write 

(*dP)(x) = t Z^(x)P i (x), 
where Pt(x) is defined by 

dP(x) = da; i P i (x). (3.6) 

The pull-back of *dP corresponding to the mapping A — > x(A) of [0, 7r] into 
M.1 n is, therefore, given by 

dz{\) := c/Ax i (A)yf(x(A))P i (x(A)) (3.7) 

for each chart whose domain contains an interval of A4q, where x l (X) and 
x 2 (A) denote the coordinates of the point x(A) in the interval. Integration 
of dz(X) now yields z(X) up to an additive constant z(0) whose value can be 
chosen freely to fit some convention. On the other hand, 

p(A) := P(x(A)). (3.8) 

The pairs (x(A),u;(A)) and (0(A), /(A)) are treated similarly to how we have 
here treated the pair (z(X), p(X)). By specifying all six of these objects, we 
are, in effect, specifying both the tangential and normal derivatives of p(x), 
/(x) and c<j(x) on the zero pressure surface. Moreover, nothing is assumed 
concerning behavior of the solution at large distances from the fluid source. 
What we are solving then is neither a Dirichlet nor a Neumann problem. 

B A Linear System for the Ernst Equation 

In a recent paper |J we described, and showed the relationships among, three 
linear systems for the vacuum Ernst equation, one inferred by the present 
authors from the Kinnersley-Chitre formalism, one developed and used ex- 
tensively by Neugebauer and a new one that the present authors have found 
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(3.5) 



valuable in the course of developing formal proofs, especially in connection 
with the hyperbolic Ernst equation.^ It is the last named one that we plan 
also to use here. It will be expressed in the form 

dQ(x,r) = A(x,r)Q(x,r) (3.9) 

for all x G M. 2 ex and r G C — K(x), where it is understood that d does not 
operate on the complex spectral parameter r (i.e., dr = 0) and where 

f r-Z(x) + P(x)* ^ f Idf(x)-Jd X (x) \ 

A(x ' r) := "I — — ){ — m^) — r 

~ (3-10) 



2/(x)' 



where 



/-(;;), *-(;_■)■ 

Ai(x, r) := [( T - Z(x)) 2 + P(x) 2 ] 1/2 , lim f^l = 1, (3.12) 

r^oo T 

and the cut lf(x) in the complex r-plane is a simple C 1 arc whose end points 
are the branch points Z(x) ±iP(x) of /i(x, r). Moreover, -K"(x) is symmetric 
with respect to the real axis, is a subset olQ 

T, ex := {z ±ip : (z, p) eVex} (3.13) 

and intercepts the real axis at a point on the same side of the closed contour 

S := {z(X) ± ip{\) : < A < tt} (3.14) 

as the point z(0) that represents the north pole n, i.e., the point of intercep- 
tion is in 

:= {zeR 1 : (z, 0) G V ex and z > z(0)}. (3.15) 
The solution of Eq. (|3.9| ) will be made unique by specifying that 



Q(n,r) := e~°^ ( J X ^ ) (3.16) 



9 In the present paper our A(x, r) and Q(x, r) are those fields that were denoted by 
rffjs(x, r) a nd (x, t) , respectively, in the cited reference. 



J See Eq. ( |2.31| ). 
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for all t G C. We reserve the option of scaling x 4 = t and choosing the 
arbitrary additive real constant in x so that £(n) = 1, whereupon <2(n, r) = 
I. To grasp the motivation behind the selection ( |3.16| ), note that the Ernst 
equation ( |2.32| ) implies that d£(z, p)/dp vanishes at p = 0. Therefore, if 
one expresses all differential forms in Eqs. fl3.9|) and (|3.10|) as functions of 
(z,p,r), 



A(.i.Q.T) = -<T.,<M~.0)-[ g J )e- 2 ^d X (z,0). 



(3.17) 



Therefore, the solution on {(z, 0) G £> e:r : z > z(0)} of Eqs. (|3.9|) and (|3.16|) 
is 



Q(z,0,r)=e 



-<T31p(z,0) 



i -x(^o) 

1 



for all z G Z ax+ and r G C. (3.18) 



Three key properties of <2(x, r) are easily deducible from Eqs. (|3.9|) and 
( |3.16| ). First, since Q(n, r) is real and A(x, r) := A(x, r*)* = A(x, r), the 
reality condition 

Q*(x,r):=[Q(x,r*)]* = Q(x.r) (3.19) 
holds. Second, since tr A(x, r) = and det Q(n, r) = 1, 



det Q(x,r) = 1. 



(3.20) 



Third, for fixed x G Ail x , Q(x, r) is a holomorphic function of r throughout 
C-K(x). 

Employing Eqs. (^T2|) , (pT3|) , (|2~29p , ( ^30|) and (plOD , one can express 
A(x, r) in the following form that no longer contains the duality operator: 



Afx,r) 



r - Z(x) //d/(x) - Jrfx(x) 



/i(x,r) 
1 



2/(x) 



0"3 



/x(x,r) 

G?x(x) 



7 ( #(x) + ^(x)dx(x)) - i J/(x)dw(x) 



J 2/(x) 
Thus, we conclude that 



0"3 

(3.21) 



0g(x(A),r) ( r-*(A) / /(A) 7 *(A)\ 
dA 1 Mx(A),r) I 2/(A) 2/(A) J 3 
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p(x(A),r) 



/ 0(A) + -w(A)*(A) 



- i J/(A>(A) 



0"3 



-J 



*(A) 



2/(A) 
and r G C - K(x{X)) 



Q(x(A),r) for all < A < vr 



(3.22) 



where Q(x(0), r) = Q(n, r) is given by Eq. (|3.16|) , where dots denote deriva- 
tives with respect to A, and where K(x(X)) lies on Eq and is 



AT(x(A)) = {z(A') ± ip(A') : < A' < A}. 



(3.23) 



We note that only the matching data (|3.2j ) and (|3.3| ) are needed in order to 
be able to write out this ordinary differential equation. We shall employ the 
convenient abbreviation 



Q(A,t) := Q(x(A),r). 



(3.24) 



Thus, both the closed contour S and the function Q(X, r) are uniquely 
determined by the matching data z(X), p(A), x(A), w(A), 0(A) and /(A), the 
differential equation ( p. 22 ) and the initial condition 



So := Q(0,r) :=e 



(73^(0) 



-x(o) 
1 



(3.25) 



C Determination of the Il-Matrix 

It is from the function Q(X, r) that one computes the monodromy matrix 
n(r) provided that the latter exists. A class of spinning fluid spheroids for 
which monodromy matrices exist will now be defined. Specifically, we shall 
henceforth consider the class of all SAPF solutions {■M.1 n ,g^ ) for which the 
following two sets of conditions hold: 

(i) The parts of the premises Q2.19| ) and ( |2.26| ) that concern M.1 n and M.f n , 
respectively, are satisfied. Also, i(A) 2 + p(A) 2 > for all A G [0, 7r]. 

(ii) As regards the function Q, there exist 2x2 matrix functions Q\ and 
0,2 with a common domain [0, 7r] x S such that EcC and is an open 
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covering of S that is either topologically equivalent to an annulus or 
is a simply connected open neighborhood of oo, E* = £, 

Q(A, r) = Qi(A, r) + /i(x(A), r)Q 2 (X, r) 

for all A G [0, tt] and r G S - # (x(A)) 1 ' ' 

and, for each r G £, Qi(A, r) and <22(A,r) are C 1 functions of A 
throughout [0,7r]. Also, for each A G [0, 7r], Qi(A, r) and dQi(X,r)/dX 
(i = 1, 2) are holomorphic functions of r throughout S. 

That completes the two sets of conditions. Conditions (ii) are required 
for the existence of the monodromy matrix. To illustrate how the conditions 
(ii) are employed in proofs, consider the fact that <2*(A,r) = <2(A,r) for all 
A G [0, tt] and r G C - K(x(A)).Q Therefore, from Eq. ( gg§ , 

[ej(A, r) - Qi(A, r)] + /*(x(A, r) [Q$(A, r) - Q 2 (A, r)] = 

for all A G [0,tt] and t G £ - K(x(X)). { ' 

For each < A < tt, analytic continuation on any simple closed orbit that 
lies in £ and that encloses one and only one of the branch points z(X)±ip(X) 
of //(x(A),r) exists and induces the replacements Qi(X,r) — ► Qj(A,r) and 
/i(x(A),t) — > — yu(x(A),r) in Eq. ( |3.27| ), whereupon one obtains Q*(A, r) = 
Qi(A,r) for all z G {1, 2}, < A < 7T and r G S - X(x(A)). The continuity 
properties of Qj then yields 

Q*(A, r) = Qi(X, r) for all i G {1, 2}, A G [0, tt] and r G E. (3.28) 

For each < A < tt, analytic continuation of ( |3.26| ) on any simple closed 
orbit in £ that encloses exactly one of the points z(X) ±ip(X) is seen to exist 
and yields 

Q'(A,r) := Q 1 (A,r)-/i(x(A),r)Q 2 (A,r). (3.29) 

Then, by letting Q^(0, r) and Q'^tt, t) denote the limits of Q'(X, r) as A — > 
and A — > 7r, respectively, we extend Eq. ( p.29| ) to all A G [0, 7f] and r G 
K(x(A)). SinceQ det Q(A,r) = 1, it follows by employing analytic contin- 
uation and the continuity properties of Qi in the manner illustrated by the 
proof of Eq. ( ggg ) that 

det Q'(X, t) = det Q(X, r) = 1 for all A G [0, tt] and r G S - i?(x(A)). 

(3.30) 



11 The deri vatio n of this fact is similar to that of Eq. ( 3.19 ). 

12 See Eq. (|3~20| ). 



14 



Furthermore, by employing analytic continuation together with the con- 
tinuities as functions of A of Qi(X,r) and dQi(X,r)/dX, one proves that 
—JQ'(X,t)J satisfies the same differential equation ( |3.22| ) as Q(A,r) for all 
A G [0, 7r] and r G £ — K(jc(X)). Thus, n(r) exists for each r G £ such that 

- JQ'(A,r)J= Q(A,r)n(r) (3.31) 

for all A G [0, 7r] and r G S — if(x(A)). By then employing Eq. ( |3.30|) and 
the fact that M T JM = (det M)J for any 2x2 matrix M, one obtains 

n(r) = -JQ(X,r) T Q'(X,r)J. (3.32) 

Therefore, with the aid of Eqs. fl3.26| ) and ( |3.29| ), we obtain the convenient 
formula 

n(r) = -jq?q'(o,t)j = -jqI^q^t) - q ], (3.33) 

from which one sees that IT(t) is holomorphic throughout £. 
From Eqs. ( |3.30| ) and ( |3.31| ) the determinantal property 

detn(r) = l (3.34) 

follows, while from Eqs. ( |3.25| ), ( |3.28| ) and ( |3.33| ) the reality property 

n*(r) = n(r) (3.35) 

follows. Moreover, from Eq. ( |3.32| ) and the fact that analytic continuation 
on a simple closed orbit that encloses one and only one of the branch points 
z{X) ±ip(X) induces Q(A,r) Q'(X,t) and Q'(X,t) Q(A,r) the symme- 
try property 

n(r) T = n(r) (3.36) 

follows. 

Upon letting A = and r = z(0) in Eqs. ( |3.26| ) and then using Eq. ( p.25| ), 
we obtain Qi(0, z(0)) = Qo Therefore, from Eq. ( [3.33j ), 

n(^(o)) = -jg%Q J. (3.37) 

It is clear from Eqs. (|3~3l ) to (gT3^) that 



n(r) is unimodular, real, symmetric and 
positive definite for all r in the maximal 
real subinterval of the real axis in S 
that contains z(0). 



(3.38) 
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D Use of II(r) to Assess Asymptotic Flatness 



It is unnecessary to construct the exterior vacuum solution in order to de- 
termine if that solution is asymptotically flat, for this can be determined by 
investigating instead the analytic properties of n(r). To help describe this 
investigation, we shall need the following concepts: 



?->EXT 
^EXT 

Sext 
Dext 



that open unbounded subset of C 
whose boundary is S , 

^EXT U S , 

{(z,p) G R 2 : (z + ip) G Z EXT }, 
{(z,p) G Sext ■ P > 0} 



and 



S ex : = ^EXT H £, 

Sex ■= {(z,p) G Sext ■ z + ip eSj, 
Vex ■= {(z,p) G S ex : p > 0}. 

There are four distinct cases: 



(3.39) 
(3.40) 
(3.41) 
(3.42) 

(3.43) 
(3.44) 
(3.45) 



(i) Suppose that S ex = Y^ext — {oo}, but 1144(7-) and r~ 1 Il34(r) have 
holomorphic extensions that cover r = 00. The transformation K 4 — > 
K 4 + A;K 3 , where k is a real number, induces n 44 (r) — > n 44 (r) and 
n 34 (r) — > n 34 (r) — 2kr. So, this transformation can be used to make 
n(r) holomorphic at r = 00, whereupon Yj ex = T>ext- Then we can 
and we do scale K 4 and select x(0) so that 



n(oo) = I. 

If it is also true that the matrix elements U a b(r) satisfy 
rn34(r)/n 4 4(r) — > as r — > 00, 



(3.46) 



(3.47) 



then the Ernst potential S of the SAV metric expressed as a function of 
Weyl canonical coordinates can be analytically extended to a domain 
that covers Sext such that £(z,—p) = £(z,p), and the restriction of 
£ to Vext yields a SAV metric without singularities that is asymptot- 
ically flat and satisfies all of the requisite matching conditions at the 
zero pressure surface. So, in this case, we obtain an asymptotically flat 
global solution. 
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(ii) Suppose that E ex = T*ext as in the preceding case, but the condition 
( |3.47|) does not hold. Then, again, £ has an analytic continuation to a 
domain that covers Sext such that £{z. — p) — £(z, p), and the restric- 
tion of £ to T>ext yields a singularity free SAV metric that satisfies all 
of the requisite matching conditions. Again, we obtain a global solution 
that some relativists would regard as asymptotically flat. However, in 
this case, in a neighborhood of spatial infinity, though 

the limit v NVT of \/ z 2 + p 2 x(z, p) as yj z 2 + p 2 — > oo (3.48) 

exists, it is not zero. For this reason, the spacetime is not asymptoti- 
cally flat in an orthodox sense. No example of this case is known, and 
it would be a shock if an example were found. 

(iii) Suppose that T> ex is a proper subspace of T>ext, and n(r) has no holo- 
morphic extension to S U S^xt (i-e., to a domain that covers T*ext), 
regardless of the choice of K4. Then £ has an analytic extension to a 
domain which covers T, ex and satisfies £(z, —p) = £(z,p). The restric- 
tion of £ to V ex yields a singularity-free SAV region that envelopes the 
fluid body and satisfies all requisite matching conditions. However, the 
full analytic continuation of £(z,p) will have at least one singularity 
on T>ext or at spatial infinity and will, therefore, not furnish a global 
solution without singularities. The Wahlquist solution may be in this 
case. 

(iv) Suppose that the holomorphic monodromy matrix n(r) does not exist 
[i.e., the conditions (ii) in Sec. |p do not hold]. Then, either our 
matching formalism is not applicable to the given SAPF or there exists 
no SAV envelope that is free of singularities (including cusps) and that 
matches the given SAPF. When our matching formalism is applicable, 
the matching SAV envelope of the given SAPF is unique (in the sense 
that any two matching SAPF envelopes will have the same full analytic 
continuation). If our matching formalism is not applicable, then there 
may exist two or more matching SAPF envelopes with different full 
analytic continuations. Criteria that would tell us [before computing 
Q(X, t)] when our formalism is applicable remain to be discovered. 

All of the conclusions that have been given above in (i) through (iv) follow 
from previous work by the authors. || |], || 
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E Constructing Exterior Solution from H(r) 

If it is desired actually to construct the SAV solution that matches the given 
SAPF solution, one proceeds to identify a Kinnersley-Chitre transformation 
matrix v(r) such thatQ 

n(r) = v(t)v(t) t , [v(t))* = v(t*), and detv(r) = 1. (3.49) 

Clearly, the matrix v (r) is defined by Eq. ( |3.49| ) only up to a transformation 

v(t) -> v(t)B(t), (3.50) 

where B(t) is an orthogonal 2x2 matrix that satisfies detS(r) = 1 and 
B*(t) = B(t). The choice of B(t) has no effect upon the SAV solution. The 
axis values of the S potential on the right side of the spinning spheroid are 
given by 



£{z,0) 



U 44 (z) 
v 33 (z) + iv 34 (z) 
-iv 43 (z) + v 44 (z) 



for all z > z(0), (3.51) 



from which, as is well known, the SAV spacetime can be constructed by many 
modern methods that are based upon Riemann-Hilbert problems. 

For example, our HHP corresponding to v(t) can be solved in two suc- 
cessive steps |J, each of which involves well known mathematics. In the first 
step, which requires only that we compute a definite integral with a given 
integrand, Minkowski space is transformed into a Weyl static spacetime. In 
the second step, which requires that we solve an ordinary Fredholm equation 
of the second kind with a given kernel and a given inhomogeneous term, the 
Weyl static spacetime is transformed into the final SAV spacetime. 



13 Some choices of the factorization IT(t) = v(t)v(t) t may only be applicable to a 
domain X — S cllt , where £* ut = T, cut and is a union of cuts in E. Each of these cuts 
crosses the real axis and its endpoints are (r-independent) isolated branch points of v(t). 
However, these branch points lead to no spacetime singularities, since the Ernst potential 
£{z,p) that is determined from v(t) will have an analytic continuation that covers T> ex . 
In summary, there may be singularities of v(t) that are not singularities of n(r); and it is 
only the singularities of il(r) that count. 
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4 Simple Illustrative Examples 



We shall illustrate all stages of the Il-matrix method with examples. In 
this discussion we shall employ geometrical units; i.e., G = 1 and c = 1. 
Moreover, the time coordinate will be scale and the arbitrary constant in x 
selected so that Q := <2(0, r) = / in Eqs. ( p.25|) and (|3.33|) . In both cases 
the matching data p(A), /(A) and u;(A) are obtained directly from the metric, 
while simple calculations yield the matching data z(\), x(A) and 0(A). 



A The Schwarzschild Interior Solution 



The first example is the Schwarzschild interior solution [jlQ 

dr 2 



1 - 2Mr 2 /R 3 
where < r < R and[^] 



(d6 2 + sin 2 6dip 2 )-e 2 ^df 



(4.1) 



^-h-\l-2Mr 2 /R 3 ) 1 / 2 



■0(x) = In 
In this chart, the pressure is given by 

p(x 



k := (1 - 2M/R) 



1/2 



(4.2) 



M 



1 - 2Mr 2 /R i - y/1 - 2M/R \ 
(47r/3)i? 3 I 3^1 -2M/R- y/T^2Mr 2 / R* ) 



(4.3) 



which is independent of 9 and which vanishes when r = R. 

Thus, for example, in the case of the Schwarzschild interior solution, the 
fields -P(x), /(x), u;(x) and x(x) are defined everywhere, and we find 



P(x) = re^ (x) sin9, /(x) = e 2 ^ (x) , u(x) = 0, x(x) = 0, 



(4.4) 



while 



★ dP(x) = - [e^(l - 2Mr 2 /R 3 ) 1/2 + k~ 1 Mr 2 /P 3 ] r sin 0d0 

+ e^(l - 2Mr 2 / R 3 )~ 1/2 cos #dr, (4.5) 
?y(x) = -k- l M{r/Rfsin9de, (4.6) 



'Note that we have scale the time coordinate t so that "0( x ) — at r = _R. 
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where we have used the relations 



*dr = -Vr 2 -2Mr 4 / J R 3 J, *d6 = dr/^/r 2 - IMr^jR?. (4.7) 
On the zero pressure surface r = R, we obtain the matching data 



Schwarzschild Interior Solution: 

z(\) = k~\R- M) cos X, p(X) = Rsm\, 

f(X) = 1, u(X) = 0, = 0, 0(A) = A; _1 M cos A + const. 



(4.8) 



The value of the constant may be chosen at will. Since only 0(A) will be 
used, there is no need to be more specific. 

Incidentally, in the case of the Schwarzschild interior solution, one can 
combine our expressions for z(X) and p(X) in the neat formula 

z(X) +ip(X) = k^Mooakifi + iX), cosh j3 := (R - M)/M, (4.9) 
which has a holomorphic extension 

z(x)+ip(yL) = k- 1 Mcosk(l3 + i0), cosh (3 := (r - M)/M, (4.10) 

where the fields z(x) and p(x) satisfy 

dz(x) = *dp(x.) and dp(x) = — *cfz(x). (4-11) 

While it is possible to employ a harmonic chart (z, p) to describe the solution, 
within the fluid the field p(x) is not the same thing as the field -P(x), which 
is quite complicated when expressed in terms of z and p. 

In the Weyl case, where o>(x) = x(x) = 0, Eq. ( |3.22j ) reduces to 

-^x~ = MxTa)^ (T3Q(A ' r) - (412) 

Introducing \I/(A, r) such that 

Q(A, r) = e-*^ 3 , (0, r) = 0, (4.13) 

we may write 



0*(A,r) _ ( T - z(XmX) + <j>(X) 
dX p(x(A),r) 
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(4.14) 



Thus, in the case of the Schwarzschild interior solution, we have 
dV(\,r) -A^MsinA 



[( T _ z (A))2 +p (A)2]l/2' 



and hence, 



*(A,r) = In 



m^-(^)-(A)-Mx(A),r) 



The function <2(A,r) is given by Eq. Q4.13Q , and 

Q'(A,r) = e -*'( A ^3 ; 

where 

m--(^)^)+Mx(A), 



*'(A,r) =ln 



Hence 
where 



n(r) :=-JQ / (0,r)J = e 2 « (T)CT3 , 



In H — In 

2 



kr — M 



kr + M / 

Therefore, we can select the K-C transformation matrix 



v(t) 



A;- 1 
k 



( kr-M \2 n 
Kkr+Mj U 

n ( kr+M \ \ 

U I kr-M I 



By Eq. ( |3.51|) the axial values of the ^-potential are given by 

kz — M 

k 2 £(z,0) = - , kz>R-M. 

v ' kz + M 



(4.15) 



(4.16) 

(4.17) 

(4.18) 

(4.19) 
(4.20) 



(4.21) 



(4.22) 



The factors k 2 and k can be suppressed by rescaling the time coordinate t and 
the spectral parameter r. We then obtain the well known axis values of the 
exterior Schwarzschild solution. The values of £(z,p) for p > are obtained 
by solving the authors' HHP or one of the other well known methods. 
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B A Winicour Dust Solution 

A somewhat more complicated example is provided by a Winicour dust 



metric 11 



ds 2 



- b2p2 '\dz 2 + dp 2 ) + p 2 dip 2 - [dt + (bp 2 /2)dipf 



(4.23) 



for which the pressure vanishes everywhere, while the energy density is given 

e(z,p) = ^b 2 e b2 " 2 /\ (4.24) 

Because the pressure vanishes everywhere, one can select any convenient 
matching surface; for example, 



z = .RcosA, p = -RsinA, 



(4.25) 



where < A < ir. Note that the requirement ^33 > yields (bR/2) 2 < 1. 

In the case of the dust metric, Z(x), -P(x), /(x), uj(x) and x( x ) are 
defined throughout the chart, and we find 

Z(x) = z, P(x) = p, /(x) = l, u^) = -hp 2 , x(*) = b(z-R), 

(4.26) 

(4.27) 



while 



b 2 p 2 dz. 



On the selected zero pressure surface ( ^.25| ), we obtain the matching data 



Winicour Dust Metric: 

z(X) = Rcos A, p(A)=i?sinA, 

/(A) = l, u{\) = -|6i? 2 sin 2 A, 

X(A) = W?(cos A - 1), 0(A) = \b 2 R 2 (cos A - f cos 3 A - §) . 



(4.28) 



When /(A) is independent of A, it is obvious from Eqs. ( |2.13| ) and ( |3.21| ) that 
those terms in A(x, r) that are proportional to the matrix 03 all vanish.0 
We find that Eq. ( |3.22D assumes the form 



dQ(X,r) 
d\ 
Q(0,r) 



-i?sin A 



m(A, r) 



&1 + J 



Q(X,t) 



(4.29) 



15 We employ units such that c = G = 1 . 

16 For this reason it is not really necessary to evaluate </>(A). 
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where 

m (A, r) := [r 2 - 2tR cos A + R 2 ) 1/2 . (4.30) 
Consider the following differential equation in the complex plane:Q 



d( 
such that 

dom T = {(C, r) G C 2 : C 7^ 00 and r £ {0, 00}} (4.32) 

and 

F(t-R,t)=I. (4.33) 

If is the function defined by the above Eqs. ( |4.31| ) to ( |4.33|) , then it is easy 
to prove with the aid of 

d -^fl^RsluX^ (4.34) 
OA m{X,T) 

and 

m (0,r) = r- J R (4.35) 

that 

Q(A,r):=^(m(A,r),r) (4.36) 
is the solution of Eqs. (f4.29|) over the domain 

{(A, r) : < A < 7T and r G C - K(x(\)) - {0, 00}}. (4.37) 

The triad of equations (|4.31|) , (|4.32|) and (|4.33f ) is collectively equivalent 
to Eq. (|4.32|) taken together with the following integral equation: 



HC, r) = I + h - J dC + £j) r) 



(4.38) 



or, equivalently, 



Q{% r ) = I + l[' d v' (l(r) + J^) S(if, r), (4.39) 



17 For fixed r, each matrix element of T satisfies the fourth order equation □ 2 J r + /3 2 J r = 
with □ := d 2 /dC 2 - (3 2 (t 2 - C 2 ) and /3 := b/(2r). We do not know if this differential 
equation or the second order equations OJ 7 ± ifiT = have been the subject of study 
elsewhere. 
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where 



and 



7 (r) := oi + J 



R 



1] := ( — (t — R) and Q(rj, r) := JF(^ + r — _R, r). 
The solution of Eq. ( |4.39| ) is given by the infinite series 



where 

and, for all n > 0, 



oo 


6) 


n 


ra=0 








r) = 


I 


rv 
/ 


7^) 


tV' 

+ J— 


/o 


T 



g n+1 (r),r) 



By mathematical induction one proves that 



Gn(v', T )- 



where 



Gn{v, T ) = —r9n(v,T), 
n\ 



k=0 



(4.40) 
(4.41) 

(4.42) 
(4.43) 
(4.44) 

(4.45) 
(4.46) 



and the coefficients g n k(T) are to be computed from the recursion relation 

71+1 



9n+l,k( T ) 

and the conditions 



n + 1 + k 



h(r)gnk(r) + Jg n ,k-i(r)} (4.47) 



9oo( T ) = I, gnkij) = if k < and if k > n. 
Thus, 0o = 



and 



0i(»7,t) = 7(t) + 



^(r,, r) = 7 (r) 2 + ( J T (r) + ± 7 (r) j) g - J^, 



(4.48) 
(4.49) 
(4.50) 
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where 



and, since Jo\ = 03, 



7(r) 



t-R 



1 13 

Ji{t) + ip{-r)J = -cr 3 - -I 



t-R 



One further proves by mathematical induction that 

g n kij) is a holomorphic function of r throughout C — {0} 

and 

T n g n k{ T ) is a holomorphic function of r 
throughout C — {00} and has the value 
(-JR) n 5 k0 at r = 0. 

From Eqs. flOfj) and flPID , 



Q(A, r) = F{m{\, r), r) = Q{m(\ 



t - r 



(4.51) 



(4.52) 



(4.53) 



(4.54) 



(4.55) 



for all (A,r) in the domain (|Q7| ). From Eqs. JO^) , Q445|) , QQ6|) , ( ^53|) , 
(|4.54j ) and the existence of the limits 



and 



lim 



lim [m(A, r) - (r - R)] = R(l - cos A) 

r — >oo 

m(A,r) - {t-R)' 



[1 — cos A) 



(4.56) 



(4.57) 



it follows as expected that 



for each A G [0, Tr],J-(m(X, r), r) has a 
holomorphic extension to all of 
C-K(x(A)), 



(4.58) 



and the equality ( |4.55|) holds for all (A, r) in the domain 

dom Q(A, r) := {(A, r) : < A < tt and r G C - i?(x(A))}. (4.59) 
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Incidentally, the reader can use Eqs. ( |4.56| ) and ( f4.57| ) to compute 



Q(A,oc) 



cxp 



-R(l -cosA)(cri + J) 



I + bR{l - cos A) 



1 




Q(A,0) 



exp 



-R(l - cos A) J 



Since <2 := Q(0, r) = / in Eq. (|3.33| ), the Pi-matrix is given by 
n(r) = -JQ(-2{t-R),t)J. 
So, from Eqs. ( |Og ), ( ggg ) and QOBD , 

[-6r(l-J2/r) 



n(r) = -J^ 

n=0 



n! 



-^(-2(r- J R),r)J 



and 



n (-2(r - i?),r) = ^fe(r) [-2(1 - R/r)}" 



k=0 



(4.60) 
(4.61) 

(4.62) 

(4.63) 
(4.64) 



Equations flOgp and ( g3D| ) imply that ^(-2(r-i?), r) and g 2 (-2(r-R), r) 
are not zero at r = oo. Further calculations reveal that g n (—2(r — R),t) 
is not zero at r = oo for all n < 10. Granting that this continues to be 
true for all values of n, we may conclude that n(r) has an isolated essential 
singularity at r = oo, and, therefore, the SAV that matches the dust metric 
at the boundary we have been considering is not asymptotically flat. A proof 
of this conclusion, however likely it seems, remains to be found. 



5 Generalizations 

The ri-matrix method can be generalized in a number of respects. The 
stationary axisymmetric source need not be of spheroidal shape. Toroidal 
sources, or multiple spheroidal sources spinning on a common axis, would 
be interesting. The source need not even be a perfect fluid. Stationary 
axisymmetric charge and current density may be involved, with resulting 
stationary axisymmetric electromagnetic fields. In principle, all such "joining 
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problems" can be handled by an extended H-matrix approach. Once the axis 
values of the the complex potentials £ and $ have been deduced from IT(t), 
the exterior electrovac fields can be constructed by solving the electrovac 
version of the authors' HHP. In this case the calculational methods devised 



by Alekseev[p^1 and Sibgatullin[|T3[ are germaine. 

We have in this first exposition of the n-matrix method avoided the math- 
ematical complications that such generalizations would entail, hoping that 
this would make it easier for the reader to appreciate the general idea behind 
this approach. 
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